Centrality and system size dependence of (multi-strange) hyperons at 40^4 and 158A 
GeV: A comparison between a binary collision and a Boltzmann + hydrodynamic 

hybrid model. 
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We present results on the centrality and system size dependence of (multi-strange) hyperons in 
Pb+Pb collisions at 40^4 and 158A GeV from the Ultra-relativistic Quantum Molecular Dynamics 
(UrQMD-v2.3) model and a coupled Boltzmann+hydrodynamics calculation. The second approach 
is realized in a hybrid fashion based on UrQMD with an intermediate hydrodynamical evolution for 
the hot and dense stage of the collision. This implementation allows a comparison of microscopic 
transport calculations with hydrodynamic simulations to explore the transition from a system that is 
not fully equilibrated such as C+C or Si+Si collisions, to a supposedly fully equilibrated system, such 
as that created in Pb+Pb reactions. The results of our calculation are compared to measurements 
of the (anti-)hyperon yields at midrapidity (|y| < 0.5) and total multiplicities performed by the 
NA49 and NA57 Collaborations at 40^4 and 158^ GeV. Furthermore, we compared our predictions 
to the centrality dependence of A, A and H~ rapidity spectra and total multiplicities at 40^4 and 
158A GeV, where possible. 
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Heavy ion collisions offer the possibility to study nu- 
clear matter under extreme conditions, at high temper- 
atures and densities, in the laboratory. This can be 
achieved by colliding nuclei, either by shooting acceler- 
ated ions at a stationary target, such as at the Alter- 
nating Gradient Synchrotron (AGS) at the Brookhaven 
National Laboratory (BNL) or the CERN Super Pro- 
ton Synchrotron (SPS), or by head-on collisions of two 
ion beams, such as at the Relativistic Heavy Ion Col- 
lider (RHIC) at the BNL and the Large Hadron Collider 
(LHC) at CERN. To achieve a "macroscopic" volume of 
excited nuclear matter, very heavy nuclei such as lead 
(Pb) or gold (An) are used. At high temperature and 
baryon density, nuclear matter is expected to melt into 
a state of free quarks and gluons, known as the quark 
gluon plasma (QGP) [l|,[3,[|. 

The strange (s) and anti-strange (s) quarks are a well- 
suited diagnostic tool to investigate the properties of the 
matter created in heavy ion reactions, because they are 
newly produced and not contained in the colliding nuclei. 
The final state distributions of strange hadrons reflect the 
strangeness production in the hot and dense stage of the 
reaction. Strangeness is a good quantum number because 
it is conserved in strong interactions. Strange quarks 
(and therefore strange hadrons) decay via the weak inter- 
action, where typical lifetimes are on the order of 10~^° 
s and so these decays are not altered on the timescale 
of a hadronizing QGP and can be well reconstructed by 
experiments. 

In Refs. d, Q the authors suggested almost 25 years 
ago that strange particle production is enhanced in the 
QGP with respect to a hadron gas and serves as a unique 
signature of the quark gluon plasma. This enhancement 



TABLE I: The energy threshold to create strange hadrons in 
hadronic interactions. 



Production reaction 


-Bthrcshoid (MeV) 


iV + 7V^A + _R'+ + iV 


^ 700 


iV-f7V-*iV + 7V + A + A 


^ 2200 


N + N ^ N + N ^S.- + -S+ 


^ 4500 




^ 5200 



can be quantifled relative to a reaction in which a tran- 
sition to a QGP phase does not take place, such as p+p 
collisions for which the system size is very small. The en- 
hancement occurs, because in the QGP different channels 
are responsible for strangeness production, leading to a 
strong reduction of the strangeness equilibration time as 
compared to the hadronic scenario. This is supported by 
the difference in threshold energies due to the fact that 
in a deconfined state, only the strange quarks have to be 
produced, rather than strange hadrons themselves. 

Hadronic transport approaches without a QGP phase 
transition usually have problems in describing the abun- 
dances of (multi-) strange baryons. In a hadronic in- 
teraction, a high threshold energy is needed to create 
strange hadrons, which is calculated from the differences 
in masses between the initial and final state particles, as 
shown in Table [J for N -I- N ^ (anti-)hyperons. Because 
of these large threshold energies (multi-)strange baryon 
production is suppressed even in initial collisions in the 
SPS energy regime « 5-20 GeV) [6, 7]. A further 

important source for (multi-strangc)hyperon production 
are multistep processes that lead to high mass resonance 
excitations and the rescatterings of initially produced an- 
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FIG. 1: Rapidity distribution of A (top) and A (bottom) from NA49 measured for different centralities in Pb + Pb collisions 
at 4:0A GeV. The closed symbols indicate measured points, whereas the open points are reflected with respect to midrapidity. 
The solid line represents calculations with UrQMD -|- hydro and the dashed line represents calculations with UrQMD. 




FIG. 2: Rapidity distribution of H~ from NA49 [3^. [3^ measured for different centralities in Pb + Pb collisions at 40 A GeV. For 
the unmeasured rapidity spectra we present the measured yield at midrapidity (|y| < 0.5). The solid line represents calculations 
with UrQMD -I- hydro and the dashed line represents calculations with UrQMD. 



tikaons. In a hadron resonance gas it takes some time to 
equilibrate strangeness and the timescale of a heavy ion 
collision might be shorter than this equilibration time. 
In regions of high net baryon density, the production of 
strange baryons is more likely, especially As produced in 
association with a kaon because this reaction has a rather 



low threshold energy. Therefore the A yields are usually 
reproduced also in nonequilibrium hadronic cascade ap- 
proaches. 

Different approaches to enhance strangeness produc- 
tion in semihadronic scenarios without a transition to a 
deconfined state have been proposed. Two different types 
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of models explain how an enhanced production of (anti-) 
hyperons can be understood. The first model [Sj] applies a 
baryon junction exchange mechanism, whereas the other 
model uses multihadronic interactions that bring the sys- 
tem to local chemical equilibrium A similar way to 
push the system into chemical equilibrium is with multi- 
ple interactions of pions 10] or the oversaturation of pion 
number [llj . Calculations in Ref. show that reaching 
an equilibrium in a pure hadronic medium, considering 
collision rates and timescales of the hadronic fireball ex- 
pansion, is only possible if hadronic multiparticlc interac- 
tions are taken into account. Recently, counterarguments 
against this conclusion were suggested within a hadronic 
model that includes Hagedorn states [l3 ]. 

The apparent enhancement of strangeness (or more 
precisely a lifting of the suppression) can be explained 
also in statistical models by the transition from a canon- 
ical ensemble, for small systems (e.g. like in p -I- p), which 
treat charge conservation exactly, to a grand canonical 
ensemble for heavier systems like Pb -I- Pb of Au + Au. 
A variation of the statistical model assuming a fully ther- 
malized core where strangeness is produced according to 
phase-space distribution surrounded by a corona which 
can be described as a superposition of elementary p-p re- 
actions has been proposed to explain the system size and 
centrality dependence of strangeness production in heavy 
ion collisions p^ . [Tsj . The equilibrium assumption seems 
to be a necessary requirement to explain the strange par- 
ticle yields. However, none of the preceding approaches 
allowed for a consistent description of both spectra and 
multiplicities of ( mult i-) strange baryons up to now. 

In this article, we attack this question by using 
a hadronic transport approach, namely the Ultra- 
relativistic Quantum Molecular Dynamics (UrQMD) in 
version 2.3 [1^, [13, [3, and compare the results to a 
hybrid approach in which a hydrodynamic evolution is 
embedded into the intermediate stage of the model dy- 
namics [l^. By doing this, we do not aim to understand 
the dynamics of the strangeness equilibration process it- 
self, but we are able to compare strangeness production 
in a full non-equilibrium approach with a dynamical ap- 
proach involving local thermal and chemical (especially 
also for strangeness) equilibrium. The advantage of this 
hybrid approach compared to a pure hydrodynamic ap- 
proach is that it allows for the direct comparison between 
non-equilibrium and locally equilibrated dynamics and 
accounts for a realistic decoupling of the particles with 
separate chemical and kinetic freeze-outs because the late 
stage rescattering is included. 

UrQMD is a microscopic many body approach to p + 
p, p -|- A and A + A interactions at relativistic ener- 
gies and is based on the covariant propagation of color 
strings, constituent quarks and diquarks accompanied by 
mesonic and baryonic degrees of freedom. The model in- 
cludes rescattering of particles, the excitation and frag- 
mentation of color strings and the formation and decay of 
hadronic resonances. Detailed information about the lat- 
est version of the UrQMD model can be found in Ref. [l^ ■ 



In the hybrid model, UrQMD serves to calculate the 
initial state of a heavy ion collision for the hydrodynami- 
cal evolution to account for the non -eq uilibrium effects in 
the very early stage of the collision 20] . With this proce- 
dure event-by-event fluctuations are naturally included. 
Furthermore, the baryon density fluctuations in one event 
are incorporated which might lead to higher strangeness 
production as discussed in Ref. [21]. The initial condi- 
tions for calculations at different centralities and for dif- 
ferent system sizes are also provided by the UrQMD ap- 
proach without further parameter adjustment. The spec- 
tators are propagated during the whole evolution in the 
hadronic cascade, whereas the other particles are mapped 
to the hydrodynamical grid. Then, the full (3 -I- 1) di- 
mensional hydrodynamic evolution proceeds employing 
the SHASTA algorithm (SHarp And Smooth Transport 
Algorithm) [l^, [2^ . For the present comparison we re- 
strict ourselves to a hadron gas equation of state which 
includes the same degrees of freedom as in UrQMD. The 
strangeness chemical potential is adjusted in the calcu- 
lation of the equation of state in such a way that the 
net-strangeness vanishes locally. In the transport ap- 
proach the local net strangeness might take a finite value, 
but this happens mostly due to the latestage rescatter- 
ings. The initial strangeness fluctuations and the influ- 
ence on the strangeness production are further discussed 
in Ref. [21] and subject to future investigations. 

The transition from the hydrodynamic phase to the 
Boltzmann dynamics, when the system is sufficiently 
diluted, takes place when the energy density e of a 
transverse slice of cells drops below 5 times the ground 
state energy density eo. This gradual transitive pro- 
cedure provides a reasonable approximation to an iso- 
eigentime criterion and has been applied successfully in 
Refs. [2^, '255. This switching criterion corresponds to a 
T — /is-configuration in which the chemical freeze-out 
is expected, approximately, at T = 170 MeV for /is = 
0. The hydrodynamic fields are then mapped to particle 
degrees of freedom via the Cooper-Frye equation on an 
isochronous hypersurface in each individual slice. For the 
production of strange baryons the corresponding chem- 
ical potentials fiB and /zg are taken into account. The 
subsequent rescatterings and final decays are calculated 
in the UrQMD model. 

Let us stress again that the present hybrid approach 
uses only a free hadron gas equation of state and has no 
phase transition included. More details about this ap- 
proach including multiplicities and spectra can be found 
in Ref. p^ . For a general comparison to similar ap- 
proaches the reader is referred to Ref. H, [13,128, 29, 3(|. 
The great advantage of this approach is that the whole 
dynamical evolution of the heavy ion reaction is gener- 
ated, and it allows us, therefore, to investigate many ob- 
servables, such as flow 12511 . HBT radii [2J|, and trans- 
verse mass spectra (2ll . l3l| at the same time and under 
the same conditions. 

In the following we show results from UrQMD and the 
hybrid approach (UrQMD + hydro) both employing only 
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FIG. 3: Rapidity distribution of A (top) and A (bottom) from NA49 ^ 33 | measured for different centralities in Pb + Pb collisions 
at 158A GeV. The closed symbols indicate measured points, whereas the open points are reflected with respect to midrapidity. 
The solid line represents calculations with UrQMD + hydro and the dashed line represents calculations with UrQMD. 




FIG. 4; Rapidity distribution of H~ from NA49 [3^ . |33| . |35| | measured for different centralities in Pb + Pb collisions at 158^1 
GeV. For the unmeasured rapidity spectra we present the measured yield at midrapidity (|y| < 0.5). The solid line represents 
calculations with UrQMD + hydro and the dashed line represents calculations with UrQMD. 



hadronic degrees of freedom. We cornpare these results 
to the measured total particle yields '42'!| and the abun- 
dances at midrapidity from the NA49 and NA57 exper- 
iments in Pb-|-Pb collisions. Furthermore, we compare 
the model predictions to the measured A, A and S~ ra- 
pidity spectra in minimum bias Pb -I- Pb and central C -I- 
C and Si + Si collisions from the NA49 Collaboration at 



40A and 158^ GeV. In addition, we present calculations 
from UrQMD and UrQMD -I- hydro for the centrality de- 
pendence of A, A and E~ rapidity spectra at Euh = 4:0 A 
GeV and Ei^h = 158A GeV and compare these to the 
measured spectra or yields at midrapidity (|y| < 0.5). 

Figures [T] and [2] show the rapidity spectra for A, A 
and E~ at E'lab = 40^1 GeV predicted from UrQMD 
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FIG. 5: The rapidity spectra of A, A and H produced in 
central C + C (left) and Si + Si (right) collisions at 158^1 
GeV [s^, [s^l . The closed symbols indicate measured points: 
open points are reflected with respect to midrapidity. The 
solid line represents calculations with UrQMD + hydro and 
the dashed line represents calculations with UrQMD. 



(dashed lines) and UrQMD + hydro (solid lines) at dif- 
ferent centrality bins and the data are presented as full 
circles (note that the open points represent the reflected 
points at midrapidity) ^32,, ,33:] . We concentrate here on 
multiplicity distributions represented by rapidity spectra 
to explore the centrality dependence of particle yields. 
UrQMD predicts a lower multiplicity compared to the 
hybrid approach. This can be understood, because at 
the transition from the hydrodynamic evolution to the 
hadronic cascade the particles are produced according 
to thermal distributions. Particle production means in 
this context the transition from macroscopic degrees of 
freedom, such as energy density and net baryon number 
density distributions to real particle degrees of freedom 
as they are propagated in the microscopic approach. In 
this manner, it is much easier to produce especially mul- 
tistrange hyperons than via resonance and string exci- 
tation. The hadronic cross sections of (multi-)strange 
baryons are small, so they do not get reequilibrated dur- 
ing the hadronic expansion phase once they are produced: 
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FIG. 6: The points correspond to the measured total yield 
by NA49 at 40 A (left) and 158^ GeV (right) [H, [H, H, [13 ■ 
The solid line represents calculations with UrQMD -I- hydro 
and the dashed line represents calculations with UrQMD. 



for example the total number of S~ at Ei^h = 40A 
GeV decreases from 2.57 directly after the hydrodynamic 
phase to 2.35 as the final value, which is far from the pure 
transport value of 0.64. The only exception are the As 
which may get absorbed even during the late stage of the 
evolution. These effects during the hadronic rescattering 
stage are usually neglected in hydrodynamic approaches 
without a hadronic cascade after the Cooper-Frye freeze- 
out. 

The A production in the pure transport approach is 
in reasonable agreement with the experimental data be- 
cause the production threshold is not as high as for the 
other hyperons because they are likely to be produced in 
nucleon-nucleon collisions in association with a kaon (see 
Table H]). Therefore the As reflect the initial baryon stop- 
ping and the rapidity distribution of the net baryon num- 
ber. In the hybrid model calculation the A production 
is higher because strangeness is chemically equilibrated. 
We overpredict the Lambda production in the hybrid ap- 
proach by 10—20 % which is similar to the deviations of 
statistical model fits to Lambda yields, as can be seen, 
for example in Ref.jsJ. In contrast to that, the nucleons 
reflect the net baryon density distribution which is ex- 
plicitly propagated in the hydrodynamic evolution. The 
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assumption of chemical equilibrium leads to enhanced A 
and S production and is in a good agreement with the ex- 
perimental data, which are shown in Figs. [3ll4l [sl. Issl. [35| 
for the highest SPS energy. There is no qualitative differ- 
ence regarding this observation between the results from 
the two beam energies. This might lead to the conclusion 
that the dynamics for strangeness production are similar 
in the whole SPS energy range. The difference between 
UrQMD (dashed line) and UrQMD-f hydro (solid hne) is 
at ii^iab = 40A GeV small for A's. 

A complementary strategy to varying the centrality of 
the collision to explore the applicability limit of the hy- 
brid model in terms of collision volume is to perform 
calculations for smaller systems. Figure [5] represents the 
rapidity spectra of A and A for central C -t- C and Si -I- 
Si collisions at 158A GeV measured by the NA49 Collab- 
oration [H, [3^ . Even if the number of participants (see 
Table in the Appendix) in central Si + Si collisions 
is similar to peripheral Pb -I- Pb collisions (33.5—43.5 % 
most central), the geometry is completely different and 
the starting time for the hydrodynamic evolution is dif- 
ferent (dependence on the radius of the nucleus R). It 
is evident that for such small systems UrQMD (dashed 
lines) describes the data better than UrQMD -I- hydro 
(solid lines). The limited system size does not allow for 
the same degree of equilibration in C -I- C and Si-|-Si 
interactions as in the case of Pb + Pb collisions for all 
shown centralities. This indicates that up to Si -I- Si reac- 
tions, equilibration cannot be achieved and the system is 
merely a (non-)trivial convolution of elementary baryon- 
baryon (and in smaller parts meson-baryon) interactions. 
From this, we also conclude that the possibility of a phase 
transition is reduced in such small systems because the 
creation of a QGP would lead to higher strangeness pro- 
duction. 

Let us next focus on a detailed study of the central- 
ity dependence of strange baryons in Pb + Pb reactions. 
Figure [H] depicts the predictions of the centrality depen- 
dence as a function of the number of wounded nucleons 
(Nw) of the total yield of A A, and S~ at i^iab = 40^ 
GeV and Ei^^ = 158 A GeV [Hi HI, H [13] ■ One clearly 
observes that for central collisions the assumption of an 
equilibrated system leads to a suitable description of the 
strangeness production within the hybrid model calcula- 
tion. To investigate if this interpretation holds also for 
smaller phase-space volumes let us turn to the particle 
production at midrapidity. 

The centrality and system size dependence of (anti-) 
hyperon production at 40^1 and 158^ GeV at midra- 
pidity (_[y| < 0.5) is shown in Fig. 7\ rneasured by the 
NA49 [HI, [H, [H, [13 and NA57 [MIm [13 Collabora- 
tions. A discrepancy between both experiments is ob- 
served. However, for a qualitative interpretation, the 
present data are sufficient and we will therefore refrain 
from speculations on the origin of the experimental in- 
consistencies. The dashed line represents calculations 
with UrQMD and the solid line depicts predictions from 
UrQMD + hydro simulations. The pure transport ap- 
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FIG. 7: Centrality dependence of the (anti-) hyperon yields 
at midrapidity (|y| < 0.5) at AOA GeV fleftj_and 158A GeV 
(right) measured by NA49 (squares) [35,'33, 36, 37] and NA57 
(circles) [38l . [39l . l40l ]. The solid line represents calculations 
with UrQMD-|-hydro and the dashed line represents calcula- 
tions with UrQMD. 



proach describes the data in very peripheral collisions 
quite well. This indicates that very peripheral Pb+Pb 
or small systems like C + C/Si + Si are not in equilib- 
rium, as discussed earlier. Moving to more central Pb 
+ Pb collisions, it is evident that the thermal produc- 
tion of strange quarks is in much better agreement with 
the measured data than obtained in the binary collision 
prescription. Especially for niultistrange hyperons it is 
visible that the production via independent strings and 
resonance excitations is not sufficient to explain the ex- 
perimentally measured yields. From this, we conclude 
that the assumption of equilibration alone allows us to 
explain a major part of the strangeness enhancement in 
central collisions of massive nuclei which does not neces- 
sarily invoke the formation of a QGP. For peripheral reac- 
tions or small collision systems a microscopic description 
of binary interactions between hadrons is favored. 

Finally, the hyperon enhancement as a function of the 
number of wounded nucleons, (Nw) ' at 158^1 GeV mea- 
sured by NA49 jFig^a)] M, M, [13 and NA57 
[Fig. Mj^)] [H, [IS [40], is shown. The strangeness en- 
hancement is defined as the yield per participant in Pb 
+ Pb collisions normalized to the yield per participant 
in p -|- p or p -|- Be interactions: 
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The centrality dependence is essentially given by the 
volume of the matter in the hydrodynamic stage at the 
transition temperature. This results in a moderate cen- 
trality dependence of the strange particle yields per par- 
ticipant from mid-peripheral to central reactions because 
the transition temperature is basically independent from 
the centrality (apart from modifications due to the slight 
change in baryon density). However, one should note 
that for very peripheral reactions an additional suppres- 
sion might be present due for the need of a canonical sta- 
tistical description [4l[ not employed here. The assump- 
tion of a grand canonical ensemble is always an infinite 
thermal bath and particle reservoir independent of cen- 
trality. For A, the absorption during the hadronic rescat- 
tering stage is reflected in the centrality dependence. A 
clear hierarchy is visible depending on the strangeness 
content. The shape of the centrality dependence is chang- 
ing from a decreasing distribution for A to an increasing 
shape for A, S~ and -I- with an enhancement for 
n~ + H"*" of approximately 10. The absorption proba- 
bility of A increases with centrality while for the other 



(multi)strange baryons the enhancement is stronger for 
larger collision volumes. UrQMD clearly underpredicts 
the effects, in all cases. 

In this article, we presented predictions for the sys- 
tem and centrality dependence of strange (anti-)hyperons 
at 40A and 158A GeV from the Ultra-relativistic Quan- 
tum Molecular Dynamics model (UrQMD) and a hybrid 
model that incorporates a (3 -I- l)-dimensional hydro- 
dynamic evolution into the UrQMD transport approach. 
In so doing this, a full non-equilibrium hadronic trans- 
port approach is compared to a hybrid model in which 
strangeness is produced according to equilibrium distri- 
butions locally. It is important to note that no explicit 
phase transition to the QGP has been taken into account 
in the present calculation. This approach takes into ac- 
count event-by-event fluctuations and is able to generate 
different centralities and system sizes without readjusting 
parameters. 

We started by showing the comparison of the central- 
ity dependence of the rapidity spectra for A, A and S~ in 
Pb + Pb collisions at Euh = 40 A GeV and E'lab = 158A 
GeV, where it was visible that due to the thermal pro- 
duction of strangeness in the UrQMD + hydro model 
a higher hyperon yield is generated. Compared to Pb 
+ Pb collisions the system seemed not fully equilibrated 
in C + C and Si + Si collisions and therefore UrQMD 
provided a better description of the measured rapidity 
spectra. The conclusion is that strangeness equilibra- 
tion is only reached in Pb + Pb collision at least in the 
shown centrality classes. This confirms the findings of 
previous studies based on statistical analysis of strange 
particle yields within an approach that includes proper 
expansion dynamics. The production of strange (anti- 
)hyperons in smaller systems like C + C and Si + Si can 
be described with a binary scattering transport approach 
like UrQMD. 
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